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Scattering  From  Snow  Backgrounds 
at  35,  98,  and  140  GHz 


1.  INTRODUCTION 

Millimeter  waves  penetrate  the  atmosphere  well  enough  to  be  of  interest  for 
military  system  applications.  Advantages  in  size,  resolution,  and  covertness 
have  stimulated  the  development  of  a  large  array  of  experimental  hardware.  Many 
of  the  required  functions  such  as  detection,  guidance,  and  communication  must  be 
carried  out  in  close  proximity  to  terrain.  Surface  scatter,  particularly  from  old, 
granular  snow,  constitutes  an  important  interaction  mechanism  with  terrain.  Be¬ 
cause  of  the  nearly  lossless  dielectric  properties  of  ice  over  the  whole  millimeter- 
wave  range,  and  because  of  particle  sizes  of  granular  snow  in  the  Mie  region,  ef¬ 
ficient  scattering  has  been  observed  over  a  wide  range  of  grazing  angles.  Snow 
granules  may  act  like  Luneberg  lenses.  That  is,  scattering  is  highly  directive  to¬ 
ward  the  direction  of  incident  radiation. 

Past  work  on  snow  and  ice  scatter  has  been  concerned  mostly  with  longer  mi- 

1  .  2 
crowaves.  Millimeter-wave  data  at  small  grazing  angles  and  at  larger  grazing 

(Received  for  publication  17  April  1984) 
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3  4  5 

angles  *  »  have  been  reported  at  35  GHz.  Some  simultaneous  measurements  were 

C 

performed  at  35  and  95  GHz.  A  good  review  of  backscatter  and  attenuation  pro¬ 
perties  of  snow  cover  at  frequencies  of  35  GHz  and  below  is  contained  in  a  report 
by  Ulaby  et  al.5 6  7 8 

In  the  RADC  program,  we  measured  millimeter-wave  interaction  with  snow 
cover  in  the  atmospheric  window  regions  of  currently  greatest  interest--35,  98, 
and  140  GHz.  Initial  results  pertained  to  the  scattering  properties  of  frozen  gran¬ 
ular  snow.  Slabs  of  this  type  could  be  cut  from  snow  cover  and  handled  without 
breaking.  A  scatterometer/transmissometer  was  designed  and  built  to  make  back¬ 
scatter,  bistatic  scatter,  and  attenuation  measurements  on  slabs.  Data  obtained  in 

O 

the  slab  mode  were  published  in  a  companion  report. 

Slabs  could  not  be  prepared  from  newly  fallen  snow  or  brittle  snow.  In  addi¬ 
tion,  the  effect  of  meltwater  on  the  scattering  properties  could  not  be  determined 
on  slabs  that  had  to  be  inverted  from  their  natural  position  during  measurement. 
Therefore,  the  investigation  was  extended  to  in-situ  snow  cover.  The  results  of 
measurements  at  35,  98,  and  140  GHz  on  in-situ  snow  during  the  winter  of  1977-78 
are  the  subject  of  this  report. 


2.  THE  IN-SITU  BACKSCATTER  PROGRAM 


The  scatterometer  employed  for  the  in-situ  measurements  was  a  modified  ver¬ 
sion  of  the  device  used  for  the  snow-slab  measurements.  The  basic  design  had 
been  proven  reliable  during  the  slab  measurements.  Staying  with  the  rugged  and 
portable  short-range  system  facilitated  comparison  between  data  obtained  in  both 
modes.  The  limited  coverage  of  the  scatterometer  precluded  wide-area  mapping  or 
imaging  of  snow-covered  terrain.  We  considered  this  unimportant,  since  we  were 
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mainly  concerned  with  the  generic  scattering  properties  of  the  medium  as  a  func¬ 
tion  of  various  microwave  and  meteorologic  parameters. 


2.1  Description  of  Scatterometer 

The  in-situ  scatterometer  is  shown  in  Figure  1.  It  differs  from  the  snow-slab 


Figure  1.  Snow  Scatterometer,  In-Situ  Backscatter  Mode  (Positioned  for  90-deg 
Grazing  Angle  Measurements) 


version  in  that  a  sample  holder  is  absent  and  the  transmitters  are  relocated  near 
the  receiver  front  end  at  the  upper  end  of  the  boom.  The  transmitter  and  receiver 
antennas  point  in  almost  parallel  directions.  Associated  power  supplies  and  other 
equipment  on  the  sled  structure  act  as  counterweights  to  the  weight  of  the  boom. 
Rotating  the  upper  part  of  the  boom  about  the  boom  midpoint  adjusts  the  grazing 
angle  between  antenna  beams  and  snow  surface.  With  the  grazing  angle  fixed  at  a 


particular  value,  the  whole  sled  assembly  rotates  around  its  vertical  axis.  During 
such  a  scan,  the  scatter  area  on  the  snow  surface  describes  a  circular  arc  whose 
radius  is  a  function  of  the  grazing  angle.  The  turntable  on  which  the  scatterometer 
rotates  can  also  be  seen  in  Figure  1.  Equipment  configurations  for  90-,  45-,  and 
15-deg  grazing  angles  appear  in  Figures  1.  2,  and  13. 


Figure  2.  Snow  Scatterometer,  In -Situ  Backscatter  Mode  (Positioned  for  15-deg 
Grazing  Angle  Measurements) 

The  block  diagram.  Figure  3,  illustrates  scatterometer  components  and  oper¬ 
ation.  In  brief,  the  transmitters  consist  of  1-kHz  squarewave  modulated  oscilla¬ 
tors  connected  to  standard-gain  horn  antennas.  A  Gunn  oscillator  generates  the 
35. 1-GHz  signal,  and  klystrons  generate  the  98. 1-  and  140. 1-GHz  signals.  Stand¬ 
ard-gain  horn  antennas  also  serve  for  reception.  The  receiver  uses  a  single-ended 
Ka-band  waveguide  mixer.  The  Gunn  local  oscillator  is  tuned  to  obtain  fundamental 
mixing  at  35  GHz,  third-order  harmonic  mixing  at  98  GHz,  and  fourth-order  har- 
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Figure  3.  Scatterometer  Block  Diagram 


monic  mixing  at  140  GHz.  The  resulting  100-MHz  IF  signal  is  amplified  and  de¬ 
tected.  Narrowband  post -detection  amplification  at  1kHz  is  performed  before  the 
signal  is  recorded  on  the  Y-axis  of  an  XY- recorder.  A  voltage  proportional  to  the 
angular  position  of  the  turntable  is  fed  to  the  recorder  X-axis.  The  snow-slab  re- 
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port  contains  a  more  detailed  description. 

2.2  Measurement  Procedures  and  Locations 

A  total  of  240  in-situ  snow  scans  were  made  over  the  period  from  10  March  to 
20  March  1978.  A  scan  produced  data  at  fixed  frequency,  polarization,  and  grazing 
angle  while  rotating  the  scatterometer  in  azimuth  through  an  arc  of  270  deg.  Scans 
were  made  at  each  of  the  three  frequencies,  at  grazing  angles  of  90,  45,  and  15 
deg,  and  at  vertical,  horizontal,  and  crossed  polarization. 

In-situ  measurements  were  carried  out  mostly  during  a  period  of  cold  nights 
and  warm  days.  Under  these  conditions,  snow  ages  through  the  so  called  melt- 
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freeze  (MF)  process.  A  day's  measurements  began  at  least  2  hours  before  sun¬ 
rise  to  obtain  a  complete  set  of  scans  while  the  snow  was  still  completely  frozen. 
Immediately  after  that,  a  second  identical  set  was  begun.  During  this  time,  free 
water  appeared  as  a  result  of  solar  heating.  When  the  free-water  content  of  the 
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snow  was  changing  rapidly,  horizontally  polarized  measurements  were  dropped  at 
|  times  to  speed  up  the  data  collection.  In  one  instance,  measurements  were  made 

at  only  one  frequency,  98  GHz. 

As  in  the  slab  mode,  the  power  ratio  was  measured  between  the  returns  from 
snow  and  from  a  large  plane  metal  calibration  plate.  Figure  4  illustrates  how 
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Figure  4.  Image  Source  Calibration  Procedure 


the  backscatter  coefficient  <7°  per  unit  area  snow  can  be  determined.  The  power 
received  from  snow  is  calculated  from  the  standard  radar  equation 

PR  =  Pf  Or  <r°  AA£ 

(I) 

<4?rrV 

with  PR  received  power 

r  distance  from  transmitter/ 

PT  transmitted  power 

receiver  to  snow 

Gt  transmit  antenna  gain 

Ae  effective  area  of  receiving  antenna 

A  snow  area  within  3-dB  points 
of  transmitted  beam.. 
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n 

-  -.i 

■  •  .1 


To  a  good  approximation, 


n&E  ^H1-2 
4  sin  6 


(2) 


where  /9E  and  /9jj  are  the  E-  and  H-plane  3-dB  beamwidths,  and  $  is  the  grazing 
angle. 

At  perpendicular  incidence  (see  Figure  4),  image  theory  yields  for  the  power 
received  from  the  metal  calibration  plate  at  range  rp 

Pbp  =  Pr  Qt  Ae.  (3) 

4  it  (2rp)2 

Combining  Eqs.  (1),  (2),  and  (3)  and  solving  for  cr°  gives 


Thus,  o°  is  calculated  from  the  measured  quantity  Pr/Prp  and  a  number  of 
known  geometric  parameters:  the  grazing  angle,  the  transmit  antenna  beamwidths, 
and  the  range  ratio. 

At  6  =  90  deg  (Figure  1),  the  antennas  passed  over  the  metal  plate  during  the 
first  40  deg  of  their  270-deg  azimuthal  scan.  The  antenna  heights  above  the  snow 
and  the  metal  plate  were  identical  in  this  case,  r g0  =  rp  =  2. 13  m.  The  scatter 
area  moved  along  an  arc  of  length  s  =  8.  8  m  during  the  270-deg  scan.  The  signal 
from  the  metal  plate  was  checked  at  the  beginning  and  at  the  end  of  a  scan,  and  the 
chart  recorder  output  was  calibrated  in  dB  below  the  metal  plate  signal,  with  a  pre¬ 
cision  attenuator  inserted  after  the  first  IF  amplifier.  At  6  =  45  deg,  the  slant 
range  to  the  snow  was  ^5  E  3.45  m,  and  s  =  16,  8  m  was  traversed  on  the  ground. 
The  corresponding  numbers  for  9  -  15  deg  are  r15  =  9. 14m  and  s  =  43.4m.  All 
like-  and  cross-polarized  measurements  on  snow  were  referenced  to  the  perpendi¬ 
cular  return  from  the  metal  plate. 

The  test  area  was  an  unobstructed,  untraveled,  essentially  flat  snow  field  of 
120  by  190  m  at  Hanscom  Air  Force  Base,  66  m  above  sea  level.  The  site  provided 
limited  shelter  for  supporting  functions  and  equipment  storage  plus  access  to  line 
power.  The  scatterometer  was  set  up  new  each  day  on  a  pedestal  of  packed  snow 
in  view  of  adjacent  areas  of  virgin  snow.  The  turntable  and  calibration  plate  were 
carefully  leveled  (Figure  1). 


Fr  4  sin# 
PRP  ^E 


3.  GROUND  TRUTH  MEASUREMENTS 


During  the  in-situ  tests,  the  following  snow  properties  were  monitored  on  a 
daily  basis: 

(1)  Depth 

(2)  Density 

(3)  Hardness 

(4)  Temperature 

(5)  Stratigraphy 

(6)  Microstructure 

(7)  Surface  Characteristics 

(8)  Free-Water  Content 

The  majority  of  these  parameters  were  determined  with  the  Snow  Observation 
Kit  developed  by  the  US  Army  Cold  Regions  Research  and  Engineering  Laboratory 
(CRREL)  and  manufactured  by  Geotest  Instrument  Corporation,  Wheeling,  Illinois. 
The  snow  was  characterized  each  day  before  the  electromagnetic  measurements 
began.  A  trench,  as  in  Figure  5,  was  established  in  an  area  of  unperturbed  snow 


near  the  backscatter  area.  The  properties  of  a  cross-section  of  snow  at  the  end  of 
the  trench  were  recorded.  Each  day,  the  trench  was  extended  to  expose  fresh 
snow.  Except  for  temperature  and  free -water  content,  the  snow  parameters  were 
not  remeasured  during  the  daily  test  period.  Air  and  snow  temperatures  were 
taken  during  each  scan.  The  free-water  content  was  measured  at  least  once  every 
hour  and  sometimes  twice  every  hour.  A  weather  diary  was  maintained  every  day, 
regardless  of  the  test  schedule,  to  record  noteworthy  changes  in  weather  condi¬ 
tions. 

3. 1  Measurement  Techniques 

Snow  parameters  (1)  through  (7)  listed  in  Section  3  were  measured  in  the  same 
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way  as  described  inthe  snow-slab  report.  Free-water  content  was  determined  only 
in  connection  with  the  in-situ  program.  Several  different  instruments  are  currently 
in  use,  and  are,  apparently,  capable  of  measuring  this  parameter  with  sufficient 
accuracy.  We  used  a  melting-water  calorimeter.  Energy  is  conserved  in  a  closed 
thermodynamic  system  undergoing  physical  changes  such  as  melting.  The  expres¬ 
sion  for  the  initial  energy  content  Ej  (calorimeter  partially  filled  with  wet  snow)  is 

Ei  -  C^MWTW  +  C*M*Ta  +  HfM*  +  <^MdTd.  (5) 

After  hot  water  has  been  added  to  melt  the  snow,  the  final  energy  content  Ef  be¬ 
comes 

Ef  =  C^(MW  +  Mg)Te  +  HfMg  +  CpMdTe.  (6) 

The  symbols  in  Eqs.  (5)  and  (6)  represent 

Cp  composite  specific  heat  of  dewar 

Cp  specific  heat  of  ice 

Cp  specific  heat  of  water 

Hf  heat  of  fusion  of  water 

Md  mass  of  dewar 

Mg  mass  of  snow  in  the  form  of  ice 

mass  of  free  water  contained  in  snow 
Mg  sum  of  Mg  and 
Mw  mass  of  hot  water  to  melt  snow 
Td  initial  temperature  of  dewar 
Te  final  equilibrium  temperature 
T  s  initial  temperature  of  snow 
Tw  temperature  of  hot  water  used  to  melt  snow 


Equating  E^  to  Ef  and  solving  for  (the  mass  of  the  frozen  part  of  the  enow), 
one  obtains 

C*MW  <TW  -  Te)  -  C^MgTe  +  (Td  -  Te) 

m;  -  - z - z -  <7) 

Hf  -  c£t8 

The  last  term  in  the  numerator  accounts  for  heat  lost  to  the  container. 

The  quantity  CpMd  was  determined  in  the  following  manner:  The  calorimeter 
(a  4-liter  glass  dewar)  was  brought  to  equilibrium  with  the  outdoor  temperature. 

A  quantity  of  water  at  room  temperature  was  filled  in,  and  the  dewar  was  agitated 
until  it  and  the  water  reached  equilibrium  temperature.  Omitting  the  snow  terms 
in  Eqs.  (5)  and  (6)  and  solving  for  CpMd,  we  find 


CpM^fTw  ~  Te) 

Te  ~Td 


Measured  results  averaged  at  C°Md  *  93  call  1C. 

Relatively  large  quantities  of  snow  (700  g)  were  used  to  determine  its  free -wa¬ 
ter  content.  With  just  enough  hot  water  (98°  C)  to  melt  the  snow  completely  but  not 
to  raise  the  final  equilibrium  temperature  much  above  0°  C,  the  heat -loss  term  in 
Eq.  (7)  was  kept  small.  During  the  experiment,  the  wide-mouthed  dewar  was  cov¬ 
ered  with  a  2-in.  cork  plate.  No  discernible  heat  loss  by  this  route  could  be  de¬ 
tected.  No  discernible  heat  loss  resulted  from  the  process  of  pouring  hot  water 
into  the  dewar.  The  overall  accuracy  of  measured  free-water  content  is  estimated 
to  be  within  2  percent. 

3.2  Snow  Conditions  During  Testa 

In-situ  backscatter  measurements  were  made  on  the  six  days  from  10  March 
to  20  March  1978.  During  this  period,  snow  conditions  varied  substantially.  On 
10  March,  the  top  layer  consisted  of  crusty  snow  that  had  been  deposited  on 
4  March.  Under  this  was  a  wide  layer  of  loose  snow  from  the  big  storm  of  6  Feb- 
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ruary.  The  bottom  layer  was  made  up  of  ice,  the  remnant  of  a  20  January  storm. 
The  surface  was  highly  irregular  because  of  wind  erosion  (Figure  6). 

Conditions  producing  strong  MF  metamorphism  remained  through  14  March. 

In  the  late  afternoon  of  this  day,  rain  began  to  fall.  The  rain,  heavy  at  times, 
ceased  in  the  evening.  The  night  of  14-15  March  was  not  cold  enough  for  the 
wet  snow  to  refreeze.  Only  late  on  16  March  did  the  temperature  drop,  and  snow 
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Figure  6.  View  of  Snow  Surface  on  11  March  1978 


began  to  fall.  On  18  March,  tests  were  conducted  over  approximately  15  cm  of 
light,  dry,  powdery  snow  that  lay  on  top  of  older  snow  and  ice.  The  surface  was 
relatively  smooth  (Figure  7).  The  next  day,  19  March,  was  warm  and  sunny.  It 


Figure  7.  View  of  Snow  Surface  on  18  March  1978 


was  followed  by  a  cold  night  that  turned  the  top  snow  cover  into  a  crusty  layer  of 
MF  snow.  The  surface  took  on  a  potholed  appearance  (Figure  8).  A  close-up  pho¬ 
tograph  of  the  MF  grains  contained  in  this  snow  is  shown  as  Figure  9.  The  small 
background  raster  is  1  by  1  mm. 

Ground  truth  data  are  compiled  in  the  Appendix  for  every  day  on  which  scatter 
measurements  took  place.  Reference  will  be  made  to  these  data  (Tables  A1-A9) 
when  discussing  the  backscatter  results  in  the  following  sections. 


Figure  8.  View  of  Snow  Surface  on  20  March  1978 


Figure  9.  Close-up  View  of  Larger  Snow  Grains  of 
20  March  1978  Shown  on  Millimeter  Grid 


4.  DATA  ANALYSIS 


Typical  raw-data  plots  recorded  during  the  backscatter  tests  are  shown  in 
Figures  10-12  for  35.  98,  and  140  GHz,  respectively.  The  plots  were  obtained 
with  the  scatterometer  set  at  a  grazing  angle  of  90  deg  and  operating  at  vertical 
polarization.  The  abscissa  in  these  figures  represents  the  azimuthal  position  of 
the  turntable  with  zero  degrees  set  at  some  arbitrary  position  beyond  the  calibra¬ 
tion  plate.  The  ordinate  is  marked  in  signal  level  in  decibel  below  that  received 
from  the  calibration  plate.  The  calibration  signal  is  seen  on  the  left  side  of  each 
plot.  Its  "flat  top"  characteristic  confirms  that  the  antenna -footprint  width  along 
the  direction  of  scan  is  substantially  narrower  than  the  width  of  the  metal  plate. 
The  irregularity  of  the  flat  top  is  attributed  to  edge  effects  and  the  limited  smooth¬ 
ness  of  the  plate. 

At  about  40  deg  from  the  beginning  of  the  sweep,  the  antenna  beams  leave  the 
calibration  plate,  and  the  fluctuating  return  from  the  snow  is  recorded  over  the  re¬ 
maining  230  deg  of  scan  angle.  Signals  in  Figures  10-12  are  similar  in  appear- 
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ance  to  those  resulting  from  backscatter  from  snow  slabs  .  The  strong  fluctua¬ 
tions  are  caused  by  the  random  distribution  of  a  large  number  of  discrete  scat- 
terers,  the  snow  granules.  Their  individual  contributions  to  the  sum  signal,  re¬ 
ceived  by  the  scatterometer  antennas,  exhibit  random  phase  that  changes  rapidly 
as  the  antennas  sweep  across  the  snow  surface. 

The  240  plots  were  analyzed  to  determine  mean  values  of  the  backscatter 
coefficient  and  other  statistical  parameters.  The  fluctuation  rate  was  measured 
as  a  function  of  signal  frequency  and  grazing  angle,  and  compared  with  a  simple 
model  of  the  scattering  process. 

4.1  Signal  Fluctuation  Rates 

During  the  slab  measurements,  the  scatterometer  antenna  beams  pointed  at  a 
fixed  location  on  the  snow  surface  that  rotated  around  an  axis  perpendicular  to  the 
beam  direction.  In  order  to  obtain  the  approximate  number  of  fluctuations  per 
degree  of  slab  rotation,  the  phasing  between  a  stationary  scatterer  at  the  axis  of 
rotation  and  one  illuminated  at  the  3-dB  point  of  the  beam  pattern  was  analyzed. 

The  majority  of  scafterers  are  closer  to  the  axis  of  rotation  than  to  the  edge  of  the 
beam.  Also,  the  product  of  transmit  and  receive  antenna  beam  patterns  weights 
the  contributions  of  the  close-in  scatterers  more  heavily  than  the  contributions  of 
those  near  the  beam  edge.  This  was  qualitatively  accounted  for  by  analyzing  phase 
changes  between  beam  edge  and  beam  center  rather  than  between  opposite  beam 
edges.  The  latter  would  have  increased  the  rate  of  change  by  a  factor  2. 

In  the  in-situ  case,  we  apply  the  same  rationale  while  adapting  the  equations 
to  the  geometry  at  hand.  Figure  13  shows  the  scatterometer  boom  configurations 


Figure  13.  Scatterometer  Positioned  at  Three  Different  Grazing  Angles 


for  measurement  at  9  =  90,  45,  and  15  deg.  The  dashed  vertical  line  represents 
the  axis  of  rotation.  The  center  of  gravity  of  the  sled  structure  stays  near  the 
center  of  the  turntable  as  the  boom's  upper  portion  is  folded  back  to  achieve  the 
45-  and  the  15-deg  look  angles.  This  leads  to  both  a  change  in  turning  radius  of 
the  antennas  and  a  change  in  height  above  the  snow  surface.  Experimentally,  the 
procedure  is  the  most  convenient  one  to  change  grazing  angles,  given  the  mechani¬ 
cal  layout  of  the  scatterometer. 

While  transmit  and  receive  horns  are  actually  mounted  side  by  side,  for  pur¬ 
poses  of  determining  the  fluctuation  rate,  they  are  assumed  to  be  co-located.  Fig¬ 
ure  14  shows  the  scatterometer  antenna  looking  down  at  the  snow  in  the  90-deg 
orientation.  Scatterers  Sj  and  S2  are  located  at  the  edge  and  at  the  center,  respec¬ 
tively,  of  the  antenna  footprint  on  the  snow  surface.  As  the  scatterometer  antenna 
moves  from  position  1  to  position  2  (distance  d),  it  is  implied  that  the  signal  scat¬ 
tered  from  S2  undergoes  no  phase  change.  We  derive  the  phase  increase  experi¬ 
enced  by  the  signal  scattered  from  Sj.  The  following  equations  are  obvious  by  in¬ 
specting  Figure  14: 

d  =  rQ  sin  a  (9) 

!  *  =  (  2  +  d2  -  2  t  d  cos  (  0  12  +  n  12).  (10) 

For  a  small  antenna  displacement  relative  to  scatterer  range,  (  >>d,  we  substi¬ 
tute  Eq.  (9)  into  (10)  and  use  a  series  approximation  for  Eq.  (10), 
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The  rate  of  range  change  with  scatterometer  rotation  around  the  vertical  axis  then 
becomes 

d  I  j/d  a  =  rQ  cos  a  sin  (0/2).  (12) 

In  the  radar  case,  the  range  change  doubles,  resulting  in  a  number  of  signal  fluc¬ 
tuations  per  incremental  angle  d  a 

d(2  i  j)/A..  =  2rQ  cos  a  sin  (  0  /2 )/  Ada.  (13) 


With  the  assumption  of  small  angles  a  and  0,  we  can  write 


This  result  is  similar  to  the  one  derived  by  Marshall  and  Hitchfield^  for  the  condi¬ 
tion  of  statistically  independent  sampling  of  the  radar  return  from  an  array  of 
scatterers  randomly  arranged  in  space.  Decorrelation  by  their  intuitive  argument 
results  when  signals  scattered  from  opposite  edges  of  the  beam  change  from  ori¬ 
ginally  equal  and  zero  phase  to  values  of  +  it  and  -  tt  by  lateral  antenna  motion 
over  distance  d.  A  more  rigorous  solution  was  published  by  Ruina, 9  10  who  gives  a 
comparable  expression  for  the  antenna  displacement  required  to  decrease  the  cor¬ 
relation  between  signals  received  to  1/e. 

Signal  fluctuation  rates  have  been  calculated  from  Eq.  (14)  for  all  wavelengths, 
scattering  angles,  and  polarizations  used  in  the  experiment.  Turning  radii  rQ  are 
shown  in  Figure  13  for  the  three  system  orientations.  Half-power  beamwidths  of 
the  transmitting  and  receiving  antennas  are  listed  in  Table  1.  Note  that  the  i.  plane 

Table  1.  E-  and  H-Plane  Beamwidths  of  Transmitting  and  Receiv¬ 
ing  Antennas 


Frequency 

Transmitting  Antenna 
Calc.  3-dB  Beamwidth 

Receiving  Antenna 
Calc.  3-dB  Beamwidth 

(GHz) 

0E(deg) 

/?H(deg) 

£E(deg) 

£H(deg) 

35 

9.1 

6.8 

13.9 

11.8 

98 

9.3 

7.0 

12. 1 

9.9 

140 

7.2 

7.2 

8.5 

7.0 

beamwidth  determines  the  fluctuation  rate  for  horizontal  polarization.  Corres¬ 
pondingly,  the  H-plane  beamwidth  is  associated  with  the  fluctuation  rate  of  verti¬ 
cally  polarized  signals.  Averages  between  transmit  and  receive  beamwidths  were 
used  for  the  calculations.  Results  for  vertical  polarization  are  plotted  as  dashed 
lines  in  Figure  15.  Also  plotted  in  Figure  15  are  averages  of  all  vertically  pola¬ 
rized  measurements  made  over  dry  snow.  The  general  slope  of  all  curves,  in¬ 
creasing  fluctuation  rate  with  increasing  grazing  angle,  is  mainly  due  to  the  rQ- 
change.  The  agreement  between  calculated  and  measurement  values  is  reasonably 

9.  Marshall,  J.  S. ,  and  Hitschfeld,  W.  (1953)  Interpretation  of  the  Fluctuating 
Echo  from  Randomly  Distributed  Scatterers,  Part  I,  Can.  J.  Phvs. 
31^:962-94. 

10.  Ruina,  J.P.  <1963)  The  fluctuation  of  radar  ground  return  due  to  antenna 
scanning,  IEEE  Trans,  on  Ant,  and  Prop.  AP-ll(No.  6):722-723. 


Figure  15.  Backscatter  Fluctuations  (2  7T-phase  Revolutions  per  De¬ 
gree  of  Scatterometer  Rotation)  vs  Grazing  Angle  (Measured - and 

Calculated  -  -  Vertical  Polarization) 


close.  The  deviations  are  not  sufficiently  systematic  to  attribute  them,  for  exam¬ 
ple,  to  the  simple  assumption  of  phasing  between  beam-edge  and  beam-center  sig¬ 
nals.  Subjective  counting  is  likely  to  be  a  greater  source  of  error.  There  is  a 
tendency  to  undercount  the  higher  fluctuation  rates  at  the  upper  two  frequencies. 
Results  for  horizontal  polarization  were  similar  and  did  not  reflect  the  differences 
calculated  in  the  E-and  H -plane  antenna  beamwidths. 

Data  from  snow  cover  containing  meltwater  showed  a  decrease  in  the  observed 
fluctuation  rate.  A  water  film  forming  on  the  surface  of  the  snow  tends  to  smooth 
it.  This  causes  a  more  specular  return  by  increasing  the  reflected,  nonfluctuating 
signal  component  arising  near  the  location  of  S2  (Figure  14).  Scattered  signal  con¬ 
tributions  from  locations  near  Sj  suffer  additional  attenuation  when  passing  obliquely 
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through  the  wet  top  layer  of  the  snow. 


I. 2  Backseatter  Statistics 

Pen  and  ink  plots  of  scattered  signals  created  some  problems  in  the  data  re¬ 
duction  process.  At  35  GHz,  extraction  of  mean  value,  standard  deviation,  and 
probability  distribution  of  signal  levels  was  readily  accomplished.  At  98  and  140 
GHz,  however,  the  scale  of  the  X-axis  (angular  position)  was  too  compressed  in 
most  cases  to  allow  visual  resolution  of  the  signal  structure  other  than  the  deter¬ 
mination  of  the  level  of  local  maxima  and  minima.  An  examination  of  Figures  10, 

II,  and  12  proves  this.  Consequently,  the  mean,  standard  deviation,  and  proba¬ 
bility  distribution  of  signal  peaks  were  determined  at  all  frequencies. 

At  35  GHz,  where  data  could  be  analyzed  completely,  a  relationship  was  es¬ 
tablished  between  the  mean  of  the  complete  signal  and  the  mean  derived  from  the 
local  peaks  alone.  It  is  assumed  that  this  relationship  holds  at  higher  frequencies, 
where  only  local  peaks  are  resolvable.  Mean  peak  levels  are  determined  by  aver¬ 
aging  over  a  large  enough  number  n  of  local  peaks  (no  more  than  20)  so  that  the 
average  changes  by  less  than  0.  2  dB  as  n  increases. 

As  an  example,  data  gathered  at  35  GHz  on  dry  snow  have  been  analyzed  for 
a  typical  day  (10  March  1978).  Figures  16,  17,  and  18  show  the  backseatter 
coefficient  in  decibel  plotted  on  probability  paper  for  vertical  polarization  and 
grazing  angles  of  90,  45,  and  15  deg.  Presented  in  this  manner,  lognormally  dis¬ 
tributed  data  fall  on  a  straight  line  and  Rayleigh -distributed  data  fall  on  a  curved 
line  of  typical  shape.  In  the  figures,  peak  data  match  straight  lines  and  complete 
data  match  Rayleigh  curves  well.  The  Rayleigh  curves  were  fitted  to  pass  through 
the  median  value  of  the  measured  a0  .  Horizontal  and  crossed-polarization  data 
resulted  in  similar  plots. 

The  Rayleigh  distribution  of  the  complete  signal  is  typical  for  a  collection  of 
randomly  positioned  scatterers,  each  contributing  roughly  equal  power.  The  log¬ 
normal  fit  of  the  peak  backseatter  data  was  somewhat  unexpected.  However, 
Gumbel, 11  in  his  book  on  the  statistics  of  extremes,  states  that  the  distribution 
of  the  extreme  values  of  a  statistic  process  quite  often  is  of  a  lognormal  nature. 

For  Rayleigh -distributed  data,  the  relationship  between  the  mean  a0  and  the 
median  <r°0  g  of  the  backseatter  coefficient  is  (in  dB) 

<J°  =  0>  5  +  l.«.  (15) 

The  cumulative  probability  curve  does  not  change  shape  for  different  sets  of  data. 


11.  Gumbel,  E,  J.  (1958)  Statistics  of  Extremes,  Columbia  University  Press, 
New  York. 


Figure  16.  Cumulative  Probability  of  <7°  PEAK  — —  A  and  a° 

A - A(10  March  1978,  35  GHz,  Vertical  Polarization,  90- 

deg  Grazing  Angle) 


and  the  standard  deviation  S  is  a  constant  (2S  =  11.3  dB).  For  lognormal  data  (in 
dB) 


a  PEAK  = 


<7°0_  5  +  0.  115  S2. 


(16) 


In  this  case,  S  is  proportional  to  the  slope  of  the  straight  line  fitted  to  the  data,  and 
is  explicitly  given  by  one  half  of  the  distance  on  the  abscissa  between  the  cumulative 


Figure  17.  Cumulative  Probability  of  ff0pEAKA - A  and 

a°A —  -  A  (10  March  1978,  35  GHz,  Vertical  Polarization,  45-deg 
Grazing  Angle) 


probabilities  of  0.  84  and  0.  16.  Table  2  lists  mean  and  mean-peak  backscatter 
coefficients  based  on  Eqs.  (15)  and  (16)  as  well  as  the  standard  deviation  of  the 
peak  backscatter  coefficient.  The  differences  between  cr°  and  <t°  found  on 

this  particular  day  are  compared  to  the  averages  of  all  35-GHz  measurements  on 
frozen  snow  at  vertical  and  horizontal  polarization.  The  overall  results  show 
o°  pgAK  *°  exceec*  a  °  by  roughly  3  dB  at  90-deg  grazing  angle,  by  2  dB  at  45- 
deg,  and  by  5  dB  at  15-deg.  As  a  firsthand  approximation,  we  assume  that  the 
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Figure  18.  Cumulative  Probability  of  o°  peak  a A  and 

- A  (10  March  1978,  35‘GHz,  Vertical  Polarization,  15-deg 

Grazing  Angle) 


mean-peak  backscatter  coefficient  is  3  dB  higher  than  the  mean  backscatter  coef¬ 
ficient  on  average,  and  that  we  can  extrapolate  from  this  behavior  at  35  GHz  to  98 
and  140  GHz. 

In  Figures  19  and  20,  the  cumulative  probability  distribution  of  <7°  is 

shown  at  98  and  140  GHz  for  the  same  day.  Again,  a  lognormal  distribution  is  a 
good  fit  to  all  data  taken  over  dry  snow.  The  mean  values  and  standard  deviations 
for  these  cases  are  presented  in  Table  3.  It  is  apparent  from  Tables  2  and  3  that 
the  backscatter  coefficient  for  dry  snow  increases  with  frequency  and  grazing 
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PROBABILITY  P  (cr>X) 


Table  2.  Mean  and  Standard  Deviation  of  Back- 
scatter  Coefficients  for  35  GHz  at  T  <  0°  C 


Grazing 

Angle 

(deg) 

Mean  Backscatter 
Coefficient 

Standard 

Deviation 

of  ®°PEAK 
S  (dB) 

<r°  (dB) 

<7°PEAK(dB) 

90 

0.5 

4.0 

3.4 

45 

-4.3 

-2.  1 

4.3 

15 

-12.8 

-6.4 

3.5 

1 

I 


■  i 
•'j 


Figure  19.  Cumulative  Probability  of  a0  peak  dO  March  1978,  98 
GHz,  Vertical  Polarization) 
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Figure  20.  Cumulative  Probability  of  o  0  peak  March  1978,  140 
GHz,  Vertical  Polarization) 


angle.  This  dependence  will  be  looked  at  in  greater  detail  in  the  next  section. 

The  data  discussed  so  far  were  collected  at  air  temperatures  of  less  than 
0°C  with  no  free  water  present.  When  the  temperature  rose  above  0°C  and  melt¬ 
water  appeared  later  during  the  day,  the  complete  signal  statistics  at  35  GHz  could 
be  established  only  at  perpendicular  incidence.  Although  more  spread  occured  in 
the  o 0  pEApT^ata  points  for  T  >  0°C  than  was  the  case  for  dry  snow  (Figure  18), 
a  lognormal  straight  line  still  provided  a  reasonable  fit  to  the  data.  However,  a 
Rayleigh  curve  no  longer  matched  the  measured  distribution  of  the  complete  signal. 
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Table  3,  Mean  and  Standard  Deviation  of  Peak  Backscatter 
Coefficient  for  98  and  140  GHz  at  T  <  0°C 


Grazing 

Angle 

(deg) 

98  GHz 

140  GHz 

^PEAK  (dB) 

S  <dB) 

^ PEAK  <dB) 

S  (dB) 

90 

4.3 

4.0 

10.7 

2.  5 

45 

3.1 

3.1 

5.  1 

2.8 

15 

-2.5 

3.0 

2.5 

4.7 

Wet-snow  <r°  pjrAK~data  March  1978  at  the  two  lower  grazing  angles  at  35 

GHz,  and  at  all  grazing  angles  at  98  and  140  GHz,  agreed  fairly  well  with  lognor¬ 
mal  distributions,  too.  Table  4  contains  O'0  pEAK*an<^  S-values  at  the  three  fre- 

Table  4.  Mean  and  Standard  Deviation  of 
Peak  Backscatter  Coefficient  for  T  >  0°C 
at  90 -deg  Grazing  Angle 


Frequency 

(GHz) 

•peak  <dB) 

S  (dB) 

35 

0.9 

5.6 

98 

2.9 

5.0 

140 

6.3 

2.7 

quencies  and  at  perpendicular  incidence  for  above-freezing  snow.  A  drop  in  the 
mean-peak  backscatter  coefficient  for  wet  snow  is  apparent  at  all  frequencies  (com¬ 
pare  Tables  2,  3,  and  4).  Again,  this  will  be  discussed  in  greater  detail  in  the 
next  section. 

1.3  Backscatter  Coefficients 

We  consider  only  peak  backscatter  coefficients  in  the  remainder  of  this  report. 
Therefore,  we  no  longer  refer  to  the  data  as  peak  data  nor  designate  the  peak  back¬ 
scatter  coefficient  as  o°  Unless  stated  otherwise,  a0  will  stand  for  mean- 

peak  values  in  what  follows.  Vertical,  horizontal,  and  cross-polarized  backscatter 


coefficients  are  distinguished  by  subscripts  VV,  HH,  and  HV.  Here  V  and  H  repre¬ 
sent  vertical  and  horizontal,  with  the  first  letter  denoting  transmitted,  and  the  se¬ 
cond  letter  denoting  received  polarization.  Backs  carter  coefficients  will  be  pre¬ 
sented  in  three  different  formats:  (a)  dry-  and  wet-snow  data  on  individual  days 
versus  grazing  angle;  (b)  averages  of  all  data  versus  grazing  angles;  and  (c)  aver¬ 
ages  of  all  data  versus  frequency. 

A  nearly  complete  set  of  measurements  on  frozen  and  melting  snow  was  ob¬ 
tained  on  two  days,  10  and  18  March  1978.  The  snow  cover  was  quite  different  on 
these  two  days.  On  10  March,  a  crusty  top  layer  had  developed  from  snow  deposit¬ 
ed  six  days  earlier.  On  18  March,  a  12-  to  15-cm  thick  layer  of  new  light  snow 
covered  the  test  area.  The  surface  was  quite  rough  on  the  earlier  day,  as  demon¬ 
strated  in  Figure  6.  It  was  relatively  smooth  on  the  latter  day,  as  seen  in  Figure 
7.  The  air  temperature  rose  rapidly  after  sunrise  on  18  March.  Measurements 
in  the  horizontally  polarized  mode  were,  therefore,  omitted  over  frozen  snow  to 
allow  the  more  rapid  collection  of  a  data  set.  Once  the  air  temperature  exceeded 
0°C,  horizontal  polarization  was  included. 

Backscatter  coefficients  at  35  GHz  are  displayed  in  Figure  21  for  10  March 
and  in  Figure  22  for  18  March.  The  general  variation  is  the  same  for  correspond¬ 
ing  curves  in  both  figures.  For  dry  snow,  ff°yy  exceeds  OdB  at  90  deg  and  is 
lower  by  10  to  15dB  at  15-deg  grazing  angle.  The  same  trend  holds  true  for 
a°  HH  (Fi8ure  21^»  although  the  drop  is  more  rapid  as  the  grazing  angle  becomes 
smaller.  Averaging  all  available  35-GHz  data  over  dry  in-situ  snow  suggests 
that  there  may  be  some  significance  to  the  more  rapid  drop-off,  but  not  to  the  dif¬ 
ference  at  perpendicular  incidence.  Cross-polarized  backscatter  coefficients  be¬ 
have  like  the  vertically  polarized  ones,  being  between  5  and  10  dB  less  in  both  fig¬ 
ures  at  each  angle  6  .  Note  the  increasing  rate  of  a  “-drop  for  dry  snow  as  the 
grazing  angle  decreases  (  6  is  plotted  on  a  linear  scale). 

The  opposite  is  true  for  the  wet -snow  data  at  this  frequency.  Free  water  pro¬ 
duced  a  dramatic  change  in  the  angular  dependence  of  o°  ,  While  the  backscat¬ 
ter  coefficient  is  only  moderately  affected  in  the  like -polarized  case  at  0  =  90  deg, 
it  drops  at  a  much  higher  rate  between  90  and  45  deg  than  previously  seen  for  dry 
snow.  Although  the  three -point  representations  of  angular  variation  of  a  0  in 
the  in-situ  mode  are  crude,  they  do  suggest  a  specular  contribution  to  wet-snow 
scatter  at  this  frequency,  similar  to  that  found  on  ice  surfaces  of  frozen  snow 
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slabs.  This  interpretation  is  supported  by  the  substantial  drop  in  <7°  ap¬ 
proaching  20dB  at  perpendicular  incidence. 

An  indication  of  a  smoother  snow  surface  on  18  March  is  seen  in  the  cr°  vv~ 
and  O  °jjjj -levels  for  wet  snow  at  90  deg.  They  are  5dB  higher  than  those  of  10 
March.  This  is  thought  to  be  the  effect  of  reflection  from  a  smooth  waterlogged 
layer  near  the  surface  of  the  snow.  Despite  the  dissimilarity  of  the  snow  on  both 


Figure  21.  O0  vs  Grazing  Angle  (10  March  1978,  35  GHz;  Dry  Snow, 
Polarization  VV  D - D,  HHD - D,  HV  D-.  -D;  Wet  Snow,  Po¬ 
larization  W  W - W,  HHW - W,  HV  W-.  -W;  Free-Water 

Content  4.  7  Percent) 


days  (Tables  A1  and  A5),  the  results  are  surprisingly  similar.  We  assume  that, 
during  the  dry  phase,  the  effect  of  the  underlying  packed  snow  is  strongly  seen. 
That  is,  the  35-GHz  signal  scattered  from  this  type  of  snow  is  not  attenuated  much 
as  it  travels  through  the  lighter  snow  on  top.  During  the  wet  phase,  snow  para¬ 
meters  other  than  the  liquid-water  content  apparently  became  less  significant  Tor 
the  scatter  process.  Tables  A 7  and  A 9  reveal  quite  similar  percentages  of  liquid 
water  during  the  measurements  in  Figures  21  and  22. 

The  curves  in  Figure  23  are  based  on  averages  of  all  35-GHz  backscatter 
coefficients  obtained  from  in-situ  snow.  The  free -water  content,  on  which  the  wet- 
snow  data  in  this  figure  are  based,  ranges  from  a  low  of  3  percent  to  a  high  of  15 
percent  by  weight.  In  view  of  the  effect  of  the  snow /water  ratio  on  backscatter, 
this  may  look  like  an  oversimplification.  However,  the  smoothing  of  the  data 
achieved  in  the  process  appears  to  benefit  the  comparison  of  all  in-situ  data  later 
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Figure  22.  a°  vs  Grazing  Angle  (18  March  1978,  35  GHz;  Dry  Snow, 
Polarization  VV  D - D,  HH  D - D,  HV  D-.  -D;  Wet  Snow,  Po¬ 
larization  W  W - W,  HH  W - W,  HV  W-.  -W;  Free-Water 

Content  4.8  to  5.3  Percent) 


on.  No  great  change  occurs  in  the  average  35 -GHz  curves  relative  to  the  ones  al¬ 
ready  displayed  for  individual  days.  Corresponding  like -polarized  ff°are  essen¬ 
tially  equal.  A  reason  for  the  slightly  higher  <7°^  for  wet  snow  and  the  slightly 
lower  for  dry  snow  is  not  apparent. 

The  specular  behavior  is  well  pronounced  for  the  overall  data.  The  relative¬ 
ly  high  return  from  wet  snow  at  15  deg  and  crossed  polarization  on  18  March  does 
not  prevail  in  the  overall  data. 

At  98  GHz,  the  backscatter  coefficient  from  dry  and  wet  snow  was  measured 
on  four  days.  On  11  March,  having  probed  the  dry  snow  at  all  frequencies,  the 
decision  was  made  to  limit  subsequent  measurements  to  only  one  frequency,  98 
GHz,  and  to  observe  the  effect  of  changing  free  water  on  backscatter  more  contin¬ 
uously.  The  98-GHz  data  of  10  March  are  shown  in  Figure  24.  Dry-snow  values 
for  like  polarization  are  higher  by  several  dB  at  perpendicular  incidence  and  conti 
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Figure  23,  Average  of  <7°  for  All  In-Situ  Data  (35  GHz;  Dry  Snow, 
Polarization  VV  D - D,  HH  D - D,  HV  D-.  -D;  Wet  Snow,  Po¬ 
larization  VV  W — W.  HH  W - W,  HV  W-.  -W;  Free-Water 

Content  3  to  15  Percent) 


nue  to  be  higher  at  the  lower  angles  in  comparison  with  averages  at  35  GHz.  Again, 
the  crossover  is  not  considered  to  be  significant.  There  is  some  indication  that 
ff°  HV  is  larger  relative  to  <7°  vv  and  <7°^  at  this  frequency.  This  is  particu¬ 
larly  the  case  at  45  deg,  and  is  also  confirmed  by  the  98-GHz  average  curves. 

Free  water  leads  to  several  characteristic  differences  between  the  results  at 
the  lower  two  frequencies.  In  Figure  24,  we  do  not  see  the  steep  drop  in  <7°Vy 
and  (away  from  0  =  90  deg)  that  was  typical  of  a  specular  scattering  pro¬ 

cess.  The  specular  contribution,  if  it  exists,  is  hardly  noticeable  for  vertical 
polarization.  For  horizontal  polarization,  the  rate  of  <7°-drop  resembles  that  of 
dry  snow.  That  is,  it  increases  with  decreasing  grazing  angle.  At  15  deg,  <7°  yv 
of  wet  snow  exceeds  of  dry  snow.  This  fact  would  be  questionable  had  it  not 
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Figure  24.  a°  vs  Grazing  Angle  (10  March  1978,  98  GHz;  Dry  Snow, 
Polarization  VV  D - D,  HH  D - D,  HV  D-.  -D;  Wet  Snow,  Po¬ 
larization  VV  W - W,  HH  W - W,  HV  W-.  -W;  Free-Water 

Content  5  Percent) 


also  been  found  on  18  March  (Figure  25),  the  other  day  on  which  horizontally  pol¬ 
arized  measurements  were  performed. 

In  contrast  to  the  35-GHz  findings,  the  cross-polarized  component  of  the  back- 
scatter  signal  is  reduced  much  less  at  98  GHz  when  the  snow  contains  meltwater 
Changes  as  large  as  20dB  were  seen  in  Figures  21  and  22,  but  the  maximum  seen 
in  Figures  24  and  25  is  about  lOdB.  Actually,  at  15  deg  on  the  18th  of  March,  the 
reduction  in  cross-polarized  backscatter  coefficient  is  only  ldB.  The  <T°j^y- 
curves  for  wet  snow  exhibit  an  angular  dependence  at  98  GHz  similar  to  the  major¬ 
ity  of  the  other  curves  at  this  frequency.  This  is  indicative  of  rough- surface  or 
volume  scatter  and  is  in  contrast  with  the  specular  appearance  of  the  wet  snow  at 
35  GHz. 

The  curves  in  Figure  26  represent  averages  of  all  data  taken  at  98  GHz  over 
wet  and  dry  snow.  Again,  as  in  Figure  23  for  35-GHz  data,  a  wide  range  of  free- 
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Figure  25.  <T°vs  Grazing  Angle  (18  March  1978,  98  GHz;  Dry  Snow, 
Polarization  VV  D-D,  HHD - D,  HV  D-.  -D;  Wet  Snow,  Polar¬ 
ization  W  W— W,  HH  W - W,  HV  W-.  -W;  Free-Water  Content 

7  to  9  percent) 


water  content  (3-16  percent)  contributed  to  the <7°-  curves.  This  should  be  kept  in 
mind  as  well  as  the  fact  that,  for  technical  reasons,  the  substantial  quantities  of 
snow  required  could  not  be  dug  from  a  shallow  layer.  Therefore,  the  free -water 
content  is  averaged  over  the  width  from  the  surface  to  the  depth  specified  in  Tab¬ 
les  A7  to  A9.  The  depth  of  snow  removal  was  chosen  differently  at  times  in  order 
to  assess  vertical  inhomogeneities  in  the  water  concentration.  Discrepancies  be¬ 
tween  snow  sampling  depth  and  millimeter-wave  penetration  depth  maybe  signifi¬ 
cant,  particularly  at  the  higher  two  frequencies. 

Visually,  it  is  obvious  that  the  average  curves  of  a°  are  less  spread  apart 
at  98  GHz  and  at  higher  levels  than  those  at  35  GHz.  For  dry  snow,  <J,0fjy  aP" 
proaches  a0  Vy  and  <T°  jjjj  more  closely  toward  the  90-deg  grazing  angle  than 
was  found  at  the  lower  frequency.  This  is  understandable  since  the  nonspecular 
scatter  process  can  provide  a  relatively  stronger  cross-polarized  return.  The 
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Figure  26.  Average  of  <f°  for  All  In-Situ  Data  (98  GHz;  Dry  Snow, 
Polarization  VV  D — D,  HH  D~  -  -D,  HV  D-.  -D;  Wet  Snow,  Po¬ 
larization  W  W - W,  HH  W - W,  HV  W-.  -W;  Free-Water 

Content  3  to  16  Percent) 


wet -snow  averages  largely  follow  the  grazing  angle  dependence  already  described 
for  individual  days.  The  segment  of  cr°yy  between  15  and  45  deg  for  melting 
snow  still  shows  a  slope  resembling  the  corresponding  35-GHz  curve  (specular 
scatter).  The  comparatively  wide  spread  and  crossover  of  a°  Vy  and  at 

the  lower  0  remains  somewhat  disturbing.  It  does  not  show  up  at  140  GHz. 

On  11  March,  data  collection  at  98  GHz  continued  all  day,  as  mentioned  before. 
The  diurnal  variations  for  the  VV-mode  and  the  HV-mode  are  seen  in  Figures  27 
and  28,  respectively.  In  each  figure,  <r°  is  plotted  as  a  function  of  time.  Air 
temperature  and  time  of  sunrise  (SR)  and  sunset  (SS)  are  marked  on  the  abscissa. 
The  ordinates  are  scaled  for  backscatter  coefficient  and  free-water  content  (FWC). 
The  heavy  dots  on  the  lower  curve  represent  measured  percentages  of  liquid  water. 
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Figure  27.  <r°  and  Free-Water  Content  (FWC)  vs  Time  of  Day  0n- 

Situ  Data  of  11  March  1978,  98  GHz,  Vertical  Polarization) 


The  curves  of  free-water  content  are  identical  in  both  figures.  The  free-water 
content  begins  to  rise  appreciably  shortly  after  0900  and  reaches  a  maximum  in 
excess  of  15  percent  before  noon.  The  high  level  of  FWC  is  maintained  to  around 
1400.  From  that  time  on,  a  steep  drop  occurs  to  around  9  percent,  which  is 
reached  at  about  1600.  Beyond  1600  and  until  the  end  of  the  measurement  at  2300, 
the  free-water  content  diminishes  at  a  slow  rate.  At  2125,  6.7  percent  of  free 
water  was  still  found  in  the  snow. 

The  period  of  heavy  melting  in  the  morning  hours  corresponds  with  a  sub¬ 
stantial  reduction  in  cr°yy  and  at  tflree  scatter  angles.  The  milli¬ 

meter-wave  data  points  during  this  period  are  too  infrequent  to  establish  an  exact 
correspondence  with  the  free-water  curve.  They  do  show  that  the  lowest  levels 
of  like-  and  cross-polarized  backscatter  coincide  with  the  period  of  highest  free¬ 
water  content  between  noon  and  1500.  An  interesting  phenomenon  is  seen  in  the 
later  part  of  the  day:  the  free-water  content  becomes  less  from  1330  on.  It  never 
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AIR  TEMPI °C) 

Figure  28.  a°  and  Free-Water  Content  (FWC)  vs  Time  of  Day  Gn- 
Situ  Data  of  11  March  1978,  98  GHz,  Cross -Polarization) 


drops  to  levels  lower  than  6.  7  percent  during  the  remainder  of  the  measurements, 
which  is  still  relatively  high.  One  would  not  be  surprised  if  the  backscatter  did  not 

g 

change  much  during  this  time.  Other  work  has  demonstrated  a  large  decrease  in 
backscatter  coefficient  even  when  as  little  as  3  percent  water  remains  in  the  snow. 
In  contrast  to  this,  the  backscatter  coefficients,  both  like-  and  cross -polarized, 
experience  a  substantial  recovery  during  the  early  period  of  disapp'  ing  free  wa¬ 
ter.  The  relatively  high  levels  of  <r°yy  and  o°Hy  are  sustained  for  several  hours 
around  sunset  before  they  begin  to  drop  again  substantially,  without  any  apparent 
clue  from  the  free -water  curve. 

An  explanation  is  possible  for  the  signal  recovery  around  1500.  As  the  snow 
melts  near  the  surface,  meltwater  percolates  down  to  the  lower  levels  of  the  snow 
cover.  Melting  caused  by  solar  heating  increases  before  noon  so  that  the  rate  of 
production  of  meltwater  exceeds  the  rate  of  downward  percolation.  Consequently, 
the  snow  near  the  surface  contains  a  higher  percentage  of  meltwater  than  the  snow 
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in  the  lower  layers.  As  the  solar  radiation  decreases  later  in  the  afternoon,  the 
percolation  rate  at  some  point  exceeds  the  rate  of  meltwater  production.  The 
free-water  content  near  the  surface  decreases  below  the  average  free-water  con¬ 
tent  sampled  by  our  calorimetric  technique.  Judging  from  the  attenuation  coeffl- 
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cient  measured  on  dry  snow  slabs  at  98  GHz,  2  to  4  dB/cm,  which  is  low  com¬ 
pared  with  the  wet  snow  under  investigation,  the  backscatter  signal  at  this  fre¬ 
quency  may  well  have  originated  in  the  top  1  or  2  cm  of  the  medium.  Possibly, 
this  pertains  only  to  times  of  extremely  strong  melting,  like  11  March.  The 
millimeter-wave  signal  sensed  snow  dryer  than  it  was  on  average. 

Over  the  3-hour  period  after  sunset  and  beyond,  free  water  continued  to  dis¬ 
appear  from  the  snow  layer  under  measurement.  By  2300,  it  became  clear  that 
the  air  temperature  would  remain  above  0°C  throughout  the  night  and  that  it  would 
not  be  possible  to  monitor  the  refreezing  process.  We  decided  at  this  time  to 
terminate  the  experiment.  Using  hindsight,  we  realize  that  data  collection  should 
have  continued  until  later  in  the  night  to  shed  light  on  the  o°yy-  and  <x°Hy-drop 
later  during  the  evening.  The  renewed  fall-off  of  the  backscatter  coefficient  at 
both  polarizations  and  at  all  grazing  angles  implies  an  increase  in  water  concen¬ 
tration  near  the  surface.  This  cannot  easily  be  explained,  since  the  sun  is  down 
and  the  air  temperature  is  falling  (Figure  28). 

The  data  in  Figures  27  and  28  are  the  result  of  a  single  day's  measurements. 
Comparison  with  the  98 -GHz  averages  in  Figure  26  shows  that,  for  dry  snow, 

<r°  yy  and  <7°HV  agree  within  2  dB  at  all  grazing  angles.  The  same  is  true  for 
<r°yV  at  perpendicular  incidence  throughout  the  melting  phase.  At  6  =  45  and  15 
deg,  a  °yy  and  a  °HV  are  lower  than  the  average  data  by  about  10  dB  during  the 
very  wet  period  in  the  morning.  They  are  within  2  to  3  dB  from  midafternoon  on. 
The  cross-polarized  <r°  at  all  grazing  angles  show  declines  during  the  wettest 
period  to  values  about  10  to  15  dB  lower  than  the  ones  in  Figure  26.  Again,  from 
midafternoon  on,  levels  move  into  the  range  shown  by  the  average  data. 

The  relatively  normal  behavior  of  <r°yy  at  6  =  90  deg  and  the  abnormally 
low  ff°yy  and  ff°Hy  in  all  other  cases  suggest  a  specular  reflection  on  the  wet 
snow.  To  the  extent  that  the  water  in  the  top  layer  appears  like  a  plane,  lossy 
dielectric  layer  to  the  millimeter -wave  signal,  no  return  is  expected  for  like  pol¬ 
arization  at  the  lower  grazing  angles.  No  cross-polarized  return  is  expected  at 
any  grazing  angle.  Whether  or  not  this  condition  was  met  by  the  snow  shown  in 
Figure  6  before  melting  is  not  known.  One  might  also  speculate  if  the  drop  in  a° 
later  during  the  evening  is  surface -structure  related  rather  than  based  on  free¬ 
water  content.  In  any  case,  these  results  point  to  the  need  for  careful  character¬ 
ization  of  snow  parameters. 

Notwithstanding  the  uncertainties  attributed  to  surface  structure,  and  the  re¬ 
quirement  to  measure  free-water  content  in  narrow  layers  close  to  the  surface. 


all  measured  reductions  in  backs catter  coefficient  A  a0  at  98  GHz  were  plotted 
against  the  percentage  of  free  water  to  obtain  proportionality  factors.  Straight 
lines  were  fitted  to  the  data  with  the  constraint  that  they  pass  through  the  origin 
(d  a°  =  OdB  for  FWC  =  0%).  The  regression  lines  average  at  A  =  0.  5  FWC 

and  d  c 0  jjy  =  FWC  with  d  <7°  in  dB  and  FWC  in  percent. 

Finally,  we  consider  the  140-GHz  data  collected  on  10  March  (Figure  29)  and 


GRAZING  ANGLE  0  (DEGREES) 

Figure  29.  a0  vs  Grazing  Angle  (10  March  1978,  140  GHz;  Dry  Snow, 
Polarization  VV  D — D,  HH  D-  -  -D,  HV  D-.  -.  -D;  Wet  Snow,  Polari¬ 
zation  VV  W — W,  HH  W - W,  HV  W-.  -.  -W;  Free-Water  Content 

5  Percent) 

18  March  (Figure  30).  Of  all  frequencies,  the  highest  backscatter  coefficients 
were  measured  at  140  GHz,  For  instance,  in  Figure  29,  a°y\/  at  90  deg  is  in  ex¬ 

cess  of  10  dB  on  10  March.  At  the  same  time,  the  effect  of  free  water  was  the 


Figure  30.  a°  vs  Grazing  Angle  (18  March  1978,  140  GHz;  Dry  Snow, 
Polarization  W  D— D,  HV  D-.  -.  -D;  Wet  Snow,  Polarization  W  W — W, 
HH  W-  -  -W,  HV  W-.  -,  -W;  Free -Water  Content  8  to  11  Percent) 


lowest  at  this  frequency.  With  a  free -water  content  of  5  percent,  the  drop  in  a° 
ranged  from  3  to  6  dB.  Similar  to  the  situation  at  98  GHz  but  unlike  that  at  35  GHz, 
the  nature  of  the  <7 -curve  did  not  change  significantly  between  dry  and  wet  snow. 

Recall  that,  on  10  March,  crusty,  old  snow  formed  the  melting  top  layer, 
whereas  on  18  March  (Figure  30),  the  upper  portion  of  the  snow  had  been  deposited 
only  two  days  earlier  and  was  essentially  fresh.  The  140-GHz  data  on  these  two 
days  reinforce  the  trend  seen  already  for  the  corresponding  98-GHz  results.  The 
spread  of  curves  is  substantially  less  on  the  second  day.  Closer  inspection  re¬ 
veals  that  this  impression  is  mostly  due  to  the  higher  cross -polarized  backscatter 
coefficient  observed  on  the  loose,  wet  snow.  Unfortunately,  only  the  curve  for 
vertically  polarized  <r°  exists  for  dry  snow  of  18  March.  Note  the  remarkable 
fact  that  it  lies  below  the  associated  wet-snow  curves  on  this  day. 

Averages  of  140-GHz  measurements  are  plotted  in  Figure  31.  Although  these 
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BACKSCATTER  COEFFICIENT  <r°(dB) 


line  with  the  and  the  corresponding  higher-frequency  curves,  which 

places  the  15 -deg  value  in  some  doubt. 

(c)  The  cross -polarized  component  changes  drastically  at  35  GHz  as  the  snow 
becomes  wet,  less  at  98  GHz,  and  almost  not  at  all  at  140  GHz. 

(d)  The  differences  between  vertical  and  horizontal  polarization  do  not  permit 
one  to  draw  systematic  conclusions,  although  there  is  an  indication  of  a° 
dropping  more  rapidly  with  grazing  angle  than  0°yy 

For  a  better  appreciation  of  their  frequency  dependence,  a°  yy  and 
have  been  redrawn  from  Figures  23,  26,  and  31.  In  Figures  32  to  34,  the  results 


Figure  32.  Average  a  0  vs  Frequency,  Grazing  Angle  of  90  deg  (Dry 
Snow,  Polarization  VV  D — D,  HV  D-.-.-D;  Wet  Snow,  Polarization 
VVW - W,  HVW-.-.-W) 


are  shown  for  one  grazing  ang?-  each,  with  frequency  as  the  independent  variable. 
The  increase  of  all  backscatter  coefficients  with  frequency  is  obvious.  A  minor 
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140 


Figure  33.  Average  a°  vs  Frequency,  Grazing  Angle  of  45  deg  (Dry 
Snow,  Polarization  VV  D — D,  HV  D-.  -D;  Wet  Snow,  Polarization 

VV  W - W,  HV  W-.  -W  ) 

deviation  is  ff°vv  at  15  deg  between  98  and  140  GHz.  At  all  grazing  angles,  dif¬ 
ferences  between  like-  and  cross -polarization  or  between  dry  and  wet  snow  de¬ 
crease  with  increasing  frequency.  In  the  majority  of  cases,  the  rate  of  <jr°-increase 
drops  with  frequency.  Frequency  is  plotted  on  a  linear  scale. 

1.1  Comparison  With  Related  Experiments 

Probably  the  earliest  measurements  of  the  backscatter  properties  of  snow 

12 

were  conducted  by  Cosgriff,  Peake,  and  Taylor.  They  measured  the  backscatter 
ceofficients  of  different  types  of  terrain  with  and  without  snow  cover.  Measure¬ 
ments  were  made  at  10,  15.  5,  and  35  GHz,  and  at  grazing  angles  ranging  from  10 


12.  Cosgriff,  R.  L.  ,  Peake,  W.  H.  ,  and  Taylor,  R.C.  (1960)  Terrain  Scattering 
Properties  for  Sensor  System  Design  (Terrain  Handbook  II).  Engineering 
Experiment  Station  Bulletin  183,  The  Ohio  State  University. 
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Figure  34.  Average  a0  vs  Frequency,  Grazing  Angle  of  15  deg  (Dry 
Snow,  Polarization  VV  D — D,  HV  D-.  -D;  Wet  Snow,  Polarization 

VV  W - W.  HV  W  -W  ) 


to  80  deg.  They  found  that,  at  35  GHz,  dry  snow  cover  led  to  an  increase  in  the 

backscatter  coefficient  over  the  coefficient  from  bare  terrain.  Wet  snow  caused 

the  coefficient  to  decrease.  Qualitatively,  these  results  agree  with  ours.  Cosgriff 

et  al  gave  insufficient  ground  truth  data  to  allow  a  quantitative  comparison.  They 

recorded  only  the  air  temperature  and  snow  depth. 

o 

Hoekstra  and  Spanogle  determined  the  backscatter  properties  of  snow  and  ice 
at  10,  35,  and  95  GHz.  Measurements  were  made  only  over  the  narrow  range  of 
grazing  angles  between  0.  3  and  1  deg.  This  is  well  outside  our  range,  15  to  90 
deg,  and  thus  precludes  a  meaningful  comparison  of  the  data.  They  do  report  10- 
GHz  backscatter  cross-sections  showing  a  rapid  drop  when  the  snow  cover  con¬ 
tained  free  water.  This  is  in  qualitative  agreement  with  our  results. 
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Currie  and  his  coworkers  of  the  Georgia  Institute  of  Technology  conducted 
a  series  of  field  measurements  on  the  backscatter  and  attenuation  properties  of 
snow  at  Steamboat  Springs,  Colorado,  in  1976.  The  backscatter  coefficient  of  dry 


and  wet  snow  was  determined  with  pulse  radars  at  35  and  95  GHz.  We  discussed 

O 

their  attenuation  measurements  in  our  snow-slab  report.  Their  backscatter 
measurements  were  made  at  15-deg  grazing  angle.  Typical  results  are  shown  in 
Table  5.  Note  that  the  values  given  for  <f°  are  termed  peak  values.  Currie  did 


Table  5.  Typical  Values  of  <?°  Measured  at  Grazing  Angle 
of  15  deg  by  Georgia  Institute  of  Technology  6 


Snow  Type 

■  ■ 

Typical  Peak  Values  of  o° 

Parallel 

Polarization 

Cross 

Polarization 

Frozen 

Metamorphic  Snow 

35 

-8dB 

-18dB 

95 

<-10dB 

<-10dB 

Wet 

Metamorphic  Snow 

35 

-15dB 

-28dB 

not  specify  whether  these  are  the  peaks  of  the  fluctuating  backscatter  signal  as  in 
our  analysis.  He  did  report,  in  agreement  with  our  results,  that  their  data  were 
best  matched  by  a  lognormal  distribution.  Their  35-GHz  data,  except  the  like- 
polarized  data  from  wet,  metamorphic  snow,  agree  quite  well  with  ours.  Meas¬ 
urements  at  95  GHz  were  limited  by  the  sensitivity  of  their  radar.  Therefore, 
backscatter  coefficients  were  obtained  only  for  frozen  snow.  Their  conclusion 
was  that  a0  at  95  GHz  is  less  than  that  at  35  GHz.  Our  Figure  34  shows  just  the 
opposite.  Currie  did  not  state  the  orientation  at  which  the  parallel  polarization 
measurements  were  made. 

The  most  comprehensive  investigations  in  the  microwave  region  on  snow  were 
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made  by  Ulaby's  group  at  the  University  of  Kansas.  *  *  *  In  1975,  they  meas¬ 
ured  backscatter  properties  of  snow  over  the  1-  to  8-GHz  range.  *  During  these 
tests,  the  snow  depth  was  typically  12  to  15  cm.  Dry  snow  had  practically  no  ef¬ 
fect  on  the  backscatter  coefficient  compared  with  the  effect  of  bare  ground.  An 
appreciable  change  was  discovered  when  the  snow  contained  free  water.  They  con¬ 
ducted  an  extended  series  of  measurements  in  1977.^'  4  These  tests  were  performed 
in  Steamboat  Springs,  Colorado,  at  frequencies  ranging  from  1  to  18  GHz  and  at 
35  GHz.  Grazing  angles  of  20  to  90  deg  were  covered.  In  addition,  radiometric 

measurements  were  carried  out  at  10.  69,  37,  and  94  GHz. 
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Our  results  complement  those  of  Ulaby  et  al.  '  *  ’  They  can  be  compared 
at  the  one  frequency  we  have  in  common,  35  GHz.  Otherwise,  there  is  no  overlap 
in  the  data.  The  snow  at  their  test  site  differed  from  our  snow  in  a  number  of 
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ways.  The  average  density  at  Steamboat  Springs  varied  from  0.21  to  0.  30g/cm  . 

3 

Ours  ranged  from  0.28  to  0. 43g/cm  .  Grain  sizes  at  their  test  site  were  smaller. 
At  Hanscom  AFB,  snow-grain  diameters  were  never  less  than  1  mm,  and,  in  most 
cases,  they  were  larger.  Due  to  the  MF  conditions  during  our  tests,  the  particles 
were  bonded  together  to  form  quite  large  conglomerates  (Figure  9).  Metamorphism 
apparently  took  place  at  Steamboat  Springs  through  a  slow  equitemperature  pro¬ 
cess,  as  Ulaby  et  al  report  grain  sizes  near  the  surface  between  0.25  and  0.5  mm. 
Only  on  the  last  day  did  they  find  particles  at  the  surface  ranging  from  1  to  2  mm 
in  diameter.  Apparently,  snow  melt  at  Steamboat  Springs  was  less  severe  than 
snow  melt  during  our  tests.  The  peak  wetness  reported  by  Ulaby  et  al  was  4. 5  per¬ 
cent.  Wetness  was  measured  in  the  top  5 -cm  layer  of  snow.  Down  to  the  same 
depth,  we  found  a  wetness  as  high  as  14.  7  percent.  A  maximum  of  17.4  percent 
was  measured  in  the  upper  7.  5-cm  layer  on  11  March.  The  differences  in  snow 
parameters  should  be  kept  in  mind  when  comparing  results. 

Figure  8-7  of  Stiles  and  Ulaby^  shows  one  set  of  <r°-data  as  a  function  of 
grazing  angle  and  polarization  over  dry  and  wet  snow  at  35  GHz.  Wet-snow  meas¬ 
urements  on  19  February  1977  were  followed  by  dry-snow  measurements  the  next 
morning.  In  Figures  35  and  36,  we  have  replotted  them  along  with  our  <r°-aver- 
ages  at  35  GHz.  Stiles  and  Ulaby  determined  larger  o°  than  we  did.  Their  values 
exceed  ours  for  dry  snow  by  between  4  and  10  dB,  depending  on  grazing  angle  and 
polarization.  For  wet  snow,  differences  cover  about  the  same  range  if  we  consid¬ 
er  like  polarization.  Cross-polarized  backscatter  coefficients  differ  by  as  much 
as  20  dB  and  more.  Both  sets  of  data  display  the  decrease  in  o°  and  a  pronounc¬ 
ed  specular  response  with  grazing  angle  as  water  forms  in  the  snow. 

We  attribute  the  observed  <7° -difference  predominantly  to  snow  structure. 

The  deviation  between  corresponding  data  points  becomes  larger  for  lower  grazing 
angles.  This  is  qualitatively  explained  by  the  less  dense,  small-grained  Colorado 
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snow  which  follows  the  volume  scatter  model  better  than  the  Massachusetts  snow. 
Massachusetts  snow  is  a  better  surface  scatterer.  The  Colorado  snow  must  also 
have  presented  a  less  mirror  like  water  surface  during  melting  or  the  cross-pol¬ 
arized  backscatter  could  not 1  *ve  been  sustained  at  the  level  shown. 

5 

Stiles  and  Ulaby  discuss  the  backscatter  response  to  grazing  angle  and  free¬ 
water  content  over  the  whole  lower  frequency  range  they  examined.  We  repeat 
their  conclusions  (in  quotation  marks)  and  compare  them  with  our  own  findings  at 
the  higher  frequencies  as  follows: 

(a)  "Snow  wetness  has  a  minor  effect  on  <7°  at  2.  6  GHz.  The  sensitivity  to 
wetness  increases  with  frequency.  " 

Our  measurements  at  35  GHz  show  a  stronger  sensitivity  of  a0  to  snow  wet¬ 
ness  than  the  University  of  Kansas  data.  We  see  the  sensitivity  decreasing  again 
at  98  and  140  GHz. 
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Figure  36.  Comparison  of  a  0  Measured  by  University  of  Kansas  at 

35  GHz  —  -  With  Our  Average  <r°  for  Wet  Snow - Taken  from 

Figure  23 


smaller  depolarization  ratio  at  35  GHz  (-7dB  at  45  deg).  The  corresponding  val¬ 
ues  are  -5dB  at  98  GHz  and -3  dB  at  140  GHz. 

(e)  "The  scatter  mechanism  seems  to  vary.  For  wet  snow,  <7°  is  governed 
by  surface  scatter  at  all  frequencies.  For  dry  snow,  volume  scattering  is  res¬ 
ponsible  for  the  a  “-behavior  at  the  higher  frequencies  and  soil  contributions  are 
important  at  the  lower  frequencies.  " 

Our  data  gave  no  indication  of  soil  influence.  At  35  GHz,  the  angular  response 
drops  off  more  rapidly  than  predicted  by  either  rough-surface  or  volume  scatter¬ 
ing.  At  the  higher  frequencies,  the  drop  is  closer  to  that  predicted  by  these  two 
models,  but  it  is  not  possible  to  distinguish  between  them. 

(f)  "Surface  roughness  effects  are  relatively  minor  (less  than  3  dB)  when  the 
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snow  is  dry. " 

The  same  conclusion  can  be  drawn  from  our  results. 

(g)  "Surface  roughness  is  very  important  when  wetness  is  present.  " 

Our  measurements  at  90  deg  showed  a  higher  return  from  wet  snow  whose 
surface  was  relatively  smooth  (compare  data  of  10  March  and  18  March).  How¬ 
ever,  at  other  grazing  angles,  we  saw  no  effect. 

(h)  "The  dynamic  range  of  <r°  due  to  wetness  increases  with  frequency  from 
about  0  dB  at  1.  2  GHz  to  as  much  as  15  dB  at  35  GHz  (for  50  deg  of  incidence)  in 
response  to  an  mv  change  from  0%  to  greater  than  1. 3%. " 

Here  m^  is  the  free-water  content  by  weight.  Our  data  agree  at  35  GHz.  At 
this  frequency,  we  measured  at  45  deg  an  average  decrease  of  16  dB  for  a0  vv 
and  12  dB  for  <r°  jjh  due  to  snow  wetness.  The  sensitivity  to  free  water  goes 
down  at  the  higher  two  frequencies.  The  reduction  in  backscatter  coefficient  at 
98  GHz  is  6  dB  for  ^°yy  and  4  dB  for  jjjj*  At  140  GHz,  the  corresponding 
numbers  are  2  dB  and  6  dB. 

(i)  "This  dynamic  range  decreases  as  9  approaches  nadir.  " 

Our  results  agree. 

(j)  "The  time  rate  of  change  of  a  0  increases  with  increasing  frequency.  This 
phenomenon  is  related  to  the  decreasing  penetration  depth  with  increasing  frequen¬ 
cy  and  the  faster  possible  rate  of  change  of  my  in  the  thinner  layer  to  which  the 
higher  frequency  0°  values  respond.  " 

We  did  not  obtain  sufficient  diurnal  information  at  35,  98,  and  140  GHz  to 
determine  response  times.  However,  we  see  no  reason  why  this  conclusion 
should  not  apply  in  our  frequency  regime. 

(k)  "ff°  response  to  wetness  is  linear  at  low  frequencies  but  becomes  non¬ 
linear  at  the  higher  frequencies.  " 

At  35  GHz,  Stiles  and  Ulaby  assume  that  <7°  is  sensitive  to  free  water  in  a 
top  layer  much  thinner  than  5  cm,  their  standard  snow  sampling  depth.  The  obser¬ 
vation  of  a  "non-linear"  dependence  of  <7° on  wetness  is  attributed  to  water  satura¬ 
tion  in  this  sensitive  top  layer,  which  coincides  with  an  average  5 -cm  free-water 
content  as  low  as  1  percent.  In  effect,  the  rate  of  <7°decrease  is  high  until  1  per¬ 
cent  of  average  free  water  is  reached,  and  much  lower  beyond  this  level. 

Our  snow  sampling  depths,  necessitated  by  the  calorimetric  equipment,  most 
likely  were  not  matched  too  well  with  the  depth  of  millimeter-wave /snow  inter¬ 
action  either.  There  is  reason  to  assume  that  near-surface  parameters  are  even 
more  relevant  at  98  and  140  GHz  than  at  35  GHz. 

(l)  "A  hysteresis  effect  is  seen  between  <7°  and  mv  in  the  top  5  cm  layer.  The 
effect  is  related  to  the  fact  that  the  response  is  governed  by  the  complete  snow 
profile  and  the  characterization  by  the  wetness  in  the  top  5  cm  layer  has  varying 
validity. " 
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Again,  "hysteresis"  appears  to  exist  because  the  wetness  of  the  directly  in¬ 
teracting  layer  could  not  be  measured.  We  may  have  seen  a  similar  effect  during 
the  period  of  exclusive  98-GHz  measurements  on  11  March  1978.  The  measure- 

O 

ments  on  dry  snow  with  presumably  greater  penetration  indicated  that  c°  varied 
little  when  a  slab  of  12.  5-cm  thickness  was  reduced  to  2.  5cm. 

(m)  "The  correlation  coefficient  and  sensitivity  of  <7°  to  changes  in  my  be¬ 
came  more  negative  with  increasing  frequency  and  angle  of  incidence. " 

Our  35-GHz  data  suggest  a  similar  dependence  on  6  .  However,  at  98  and 
140  GHz,  it  is  not  clear  that  free-water  effects  change  the  same  way  with  angle. 
We  definitely  observe  a  frequency  dependence  above  35  GHz  reversed  from  the  one 
found  by  Stiles  and  Ulaby. 


5.  CONCLUSIONS 

Backscatter  measurements  at  35,  98,  and  140-GHz  over  undisturbed  in-situ 
snow  cover  were  conducted  to  study  the  effect  of  liquid  water  on  the  millimeter- 
wave  backscatter  coefficient,  and  to  complement  previous  results  on  frozen  snow 
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slabs.  Since  relatively  steep  depression  angles  appeared  to  be  of  interest  from 
an  applications  viewpoint,  grazing  angles  of  90,  45,  and  15  deg  were  selected. 
Comparability  of  the  slab  measurements  and  dynamic -range  limitations  in  the  in¬ 
strumentation  contributed  to  this  choice.  Although  areas  covered  by  individual 
measurements  were  substantially  larger  in  the  in-situ  case  than  in  the  slab  case, 
this  was  not  an  attempt  at  large-scale  mapping  or  imaging  of  the  snow  surface 
which  might  be  typical  for  an  airborne  radar  system.  The  measurements  were 
conducted  over  level  terrain,  and  the  snow  appeared  to  be  uniform  within  the  re¬ 
gion  covered  by  the  antenna  beams.  Backscatter  data  recorded  during  the  circular 
sweeps  of  the  3-frequency  scatterometer  across  the  snow  surface  were  therefore 
averaged  over  each  sweep.  Structurally,  the  snow  varied  from  fresh  to  metamor- 
phic,  almost  4  weeks  old.  The  total  depth  was  never  less  than  30  cm.  Based  on 
attenuation  coefficients  measured  on  slabs,  it  appears  that  scatter  from  the  under¬ 
lying  terrain  can  be  neglected  even  at  35  GHz,  where  the  deepest  penetration  oc¬ 
curred. 

During  the  test  period,  the  snow  was  usually  frozen  solid  before  sunrise. 
Within  1  or  2  hours  after  sunrise,  free  water  could  be  detected,  A  melting-water 
calorimeter  provided  the  free-water  data  reported  here.  It  must  be  understood 
that  measuring  free-water  content  by  fractional  weight  and  interpreting  its  effect 
on  backscatter  is  difficult.  As  indicated  in  the  body  of  the  report,  free-water 
generation  in  the  snow  and  vertical  transport  of  free  water  through  the  snow  is  a 
complicated  process.  The  need  to  probe  snow  consistency  within  narrow  layers  is 
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somewhat  in  conflict  with  the  requirement  for  substantial  quantities  of  snow  in  or¬ 
der  to  reduce  calorimeter  error.  Free -water  concentrations  as  measured  show 
some  nonsystematic  variation,  and,  hence,  scatter  when  correlated  with  backscat- 
ter  coefficients.  The  adequate  characterization  of  wet  snow  is  a  problem  in  other 
areas  of  research,  too.  Notwithstanding  these  limitations,  the  general  trends  ob¬ 
served  for  the  influence  of  meltwater  on  snow  backscatter  are  believed  to  be  cor¬ 
rect. 

In  summary,  the  increase  in  the  like -polarized  backscatter  coefficient  with 
frequency,  found  during  the  slab  measurements  on  dry  snow,  was  confirmed  at  all 
grazing  angles.  Values  ranged  from  a  minimum  of  -9dB  at  35  GHz  and  15-deg 
grazing  angle  to  well  over  OdB  at  all  frequencies  (8dB  to  140  GHz)  't  90-deg  graz¬ 
ing  angle.  The  cross-polarized  backscatter  coefficient,  following  the  same  gener¬ 
al  trend,  spread  from  a  low  of  -16dB  at  35  GHz  and  15  deg  to  a  high  of  2dB  at  140 
GHz  and  90  deg.  With  increasing  frequency,  volume  or  rough-surface  scatter  ap¬ 
peared  to  prevail  over  specular  reflection.  This  was  also  apparent  when  liquid 
water  formed  in  the  snow.  While  in  almost  all  cases  like -polarized  backscatter  co¬ 
efficients  were  lower  for  wet  snow  than  for  dry  snow,  the  drop  became  less  signi¬ 
ficant  with  higher  frequency.  Actually,  at  140  GHz  and  perpendicular  incidence, 
no  reduction  was  found.  Cross-polarized  backscatter  coefficients  of  wet  snow  at 
90-deg  grazing  angle  ranged  from  a  low  of  -24dB  at  35  GHz  to  a  high  of  OdB  at  140 
GHz.  At  lower  grazing  angles,  the  change  was  somewhat  less  dramatic.  Overall, 
the  frequency  dependence  for  wet  snow  was  more  pronounced  than  for  dry  snow. 


References 


1.  Ulaby,  F.T.,  Stiles,  W.  H. ,  Dellwig,  L.  F. ,  and  Hanson,  B.  C.  (1977)  Experi¬ 

ments  on  the  radar  backscatter  of  snow,  IEEE  Trans,  on  Geoscience 
Electronics.  GE-15  (No.  4):  185-189. 

2.  Hoekstra,  P. ,  and  Spanogle,  D.  (1972)  Backscatter  from  snow  and  ice  sur¬ 

faces  at  near  incident  angles,  IEEE  Trans,  on  Ant,  and  Prop. ,  AP-20 
(No.  6):  788-790. 

3.  Stiles,  W.H.,  Hanson,  B.C.,  and  Ulaby,  F.T.  (1977)  Microwave  Remote 

Sensing  of  Snow,  Experiment  Description  and  Preliminary  Results.  Univer¬ 
sity  of  Kansas  RSL  Tech.  Hep.  340-1. 

4.  Ulaby,  F.T.,  and  Stiles,  W.  H.  (1978)  The  Active  and  Passive  Microwave 

Response  to  Snow  Parameters.  Part  I:  Wetness  Part  II:  Water  Equivalent 
of  Dry  Snow.  University  of  kansas  RSL  Tech.  Rep.  340-2. 

5.  Stiles,  W.H.,  and  Ulaby,  F.T.  (1980)  Microwave  Remote  Sensing  of  Snow- 

packs.  NASA  Contr.  Rep.  3263,  Contract  NAS5- 23777. 

6.  Currie,  N.  C.  (1977)  Private  communication. 

7.  Ulaby,  F.T.,  Fung,  A.K..  and  Stiles,  W.  H.  (1978)  Backscatter  and  Emission 

of  Snow:  Literature  Review  and  Recommendations  for  Future  Investigations. 
University  of  Kansas  RSL  Rep.  369-1. 

8.  Lammers,  U.H.W.,  Hayes,  D.T.,  and  Marr,  R.A.  (1981)  Millimeter  Wave 

Scatter  and  Attenuation  Measurements  on  Snow  Slabs.  RADC-TR-81-88, 

AD  A 112595. 

9.  Marshall,  J.  S.,  and  Hitschfeld,  W.  (1953)  Interpretation  of  the  Fluctuating 

Echo  from  Randomly  Distributed  Scatterers,  Part  I,  Can.  J,  Ph.ys, 
31:962-94. 

10.  Ruina,  J.  P,  (1963)  The  fluctuation  of  radar  ground  return  due  to  antenna 

scanning,  IEEE  Trans,  on  Ant,  and  Prop.  AP-ll(No.  6):722-723. 

11.  Gumbel,  E.J.  (1958)  Statistics  of  Extremes.  Columbia  University  Press, 

New  York, 


r.* 


i 


12.  Cosgriff,  R.  L.,  Peake,  W.  H. ,  and  Taylor,  R.C.  (1960)  Terrain  Scattering 
Properties  for  Sensor  System  Design  (Terrain  Handbook  U).  Engineering 
Experiment  Station  Bulletin  181,  The  Ohio  State  University. 


56 


Appendix  A 

Snow  Data 

Physical  parameters  of  the  snow  cover  are  presented  in  this  section.  One 
combination  table /diagram  is  given  for  every  day  that  millimeter-wave  measure¬ 
ments  were  conducted  (Tables  A1  through  A6).  The  data  shown  in  these  tables 
characterize  the  snow  before  the  time  the  millimeter-wave  measurements  were 
begun.  A  standard  form  of  presentation  has  been  used.  The  vertical  scale  at  the 
center  of  the  form  gives  the  height  above  ground.  The  stratigraphy  of  the  snow 
cover  is  shown  by  horizontal  lines  drawn  between  the  left  edge  of  the  form  and  the 
height  scale.  Most  of  the  column  headings  are  self-explanatory.  The  "Layer" 
column  normally  contains  a  brief  description  of  the  layer  such  as  hard,  icy, 
slightly  wet,  etc.  In  the  "Density"  column,  a  tick  mark  is  placed  to  show  the 
height  at  which  the  measurements  were  made.  The  air  temperature,  indicated  by 
the  symbol  (J)  on  the  temperature  plot,  was  obtained  at  a  position  just  above  the 
surface  of  the  snow.  The  "Comments"  column  includes  such  information  as  sur¬ 
face  hardness,  weather  data,  data  on  the  history  of  the  snow,  and,  in  some  cases, 
additional  information  on  the  stratigraphy  of  the  snow.  The  snow  was  mostly 
frozen  when  the  pre-dawn  data  shown  in  Tables  A1  through  A6  were  collected. 

A  corresponding  set  of  tables  (A7  through  A9)  contains  the  diurnal  variation  of 
snow  and  air  temperatures  as  well  as  the  change  in  free-water  content.  These 
data  were  recorded  at  intervals  of  one-half  hour  or  greater,  as  the  change  in  con¬ 
ditions  warranted.  "Depth  of  Sample"  means  that  snow  from  the  surface  down  to 
the  indicated  level  was  used  to  determine  the  free-water  content. 
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